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We present a method for reahzing efficiently Grover's search algorithm in an array of coupled 
cavities doped with three-level atoms. We show that by encoding information in the lowest two 
ground states of the dopants and through the application of appropriately tuned global laser fields, 
the refiection operator needed for the quantum search algorithm can be realized in a single physical 
operation. Thus, the time steps in which Grover's search can be implemented become equal to the 
mathematical steps ~ 0{y/N), where A*' is the size of the register. We study the robustness of 
the implementation against errors due to photon loss and fluctuations in the cavity frequencies and 
atom-photon coupling constants. 
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The potential to execute certain types of algorithms 
much more efficiently than the corresponding classical 
counterparts is one major reason why quantum comput- 
ers were initially proposed. To date, several important 
quantum algorithms have been invented: the two most 
well known ones being the quantum search and the prime 
number factorization ^] . To realize a physical implemen- 
tation of a quantum algorithm, it is important that one 
can encode the basic units of quantum information, i.e 
the qubits, initialize them to some suitable inputs, per- 
form an adequate set of unitary operations and then fi- 
nally read the output. Many physical implementations 
have been proposed and performed based on trapped 
ions, cold atoms and solid-state systems Q with vary- 
ing degree of success for realizing few-qubit applications. 
More recently, coupled cavity arrays have been proposed 
as a new hybrid light-matter system for quantum sim- 
ulation applications. This system was initially studied 
for conditionalphoton phase gates @ and later for Mott 
transitions [^-l6|. Following the initial works, a number 
of papers appeared studying the details of the polari- 
tonic many-body state and simulations of more complex 
spin models were proposed 0-13 ■ More recently driven 
arrays were considered towards the production and co- 
herent control of steady state entanglement [l^ under 
realistic dissipation parameters. Also, an analogy with 
Josephson oscillations was shown and the many body 
properties of the driven array have been recently stud- 
ied [H]. 

In the current work we present for the first time a 
method for realizing Grover's search algorithm for an ar- 
ray of coupled cavities in an efficient way. We start by 
noting that the basic Grover iteration is executed in two 
steps - first, the oracle O is applied to mark the searched 
state by flipping its phase, and then a global reflection 
operation about the mean G{W) = 1 — 2|VF)(Vl^| is per- 
formed, where W is the A^-qubit W state. We will show 
that by exploiting the natural evolution of our system 
we can implement the reflection operation G{W) - also 
known as a quantum mirror or Householder reflection [lij 
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FIG. 1: An array of coupled cavities via hopping photons 
between neighbouring sites at a rate J doped with three-level 
atomic lambda systems. 



- in only one physical step. The quantum reflection, to- 
gether with the oracle operation O, comprise the Grover 
logical step which is performed ^ 0(2-^) times to find 
the searched entry. The current techniques for perform- 
ing these transformations utilize sequences of single- and 
two- qubit gates in order to perform an A-qubit quan- 
tum gate [20| . Hence, the number of operations scale as 
^ ©(A^) with the size of the register. Here, we will ex- 
emplify how to perform the required A^-qubit reflection 
gate for the Grover iteration in only one operational step. 
In this sense the required physical time steps are signifi- 
cantly reduced and become equivalent to the number of 
mathematical steps ~ 0{ ^^ ). 

Consider a linear chain of A^ cavities fabricated in such 
a way that the spatial profile of their cavity modes over- 
lap and photons can hop between neighbouring sites. Let 
each cavity be doped with a three-level atomic system 
with a lambda internal energy configuration, as shown in 
Fig. 1 and let us assume that the A^ logical qubits, which 
we utilize for the quantum register, are encoded in the 
ground states |0)fe and |l)fc of the individual atoms fssj. 
The |0)fc O \e)k transition of each atom is coupled to 
the respective cavity mode with coupling strength g and 
detuning A = wq — Wc being the difference between the 
Bohr transition frequency of the atom ujq and the cavity 
frequency Wc- Additional laser fields drive the |l)fe 4-> \e)k 
transitions with Rabi frequency f2 and are detuned from 
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the transition frequency again by A. The total Hamilto- 
nian which governs the evolution of the system consists 
of three parts - a free energy term, an interaction term 
describing the coupling between the atoms and the re- 
spective cavity modes, and a hopping term. Using the 
dipole and the rotating wave approximations we can ex- 
press the system Hamiltonian in an appropriately chosen 
interaction picture as 

N N 

,_,totai ^ J2^caiak + jJ2i('W+i+il.c.) 

k=l k=l 
N 

fe=i 
1 ^ 

+ -17^(e-^*|l),(eU-|-H.c.), (1) 

fe=i 

where we have adopted the convention h — 1. Without 
loss of generality we can assume that and g are real, by 
including their phases in the atomic and cavity photon 
states. 

For the one step realization of the Grover search algo- 
rithm we will exploit the single-excitation subspace and 
initialize all qubits in their ground states \0)k while pop- 
ulating the lowest energy common photonic mode with 
one photon. Hence, the available Hilbcrt subspace is 
(A^ -I- l)-dimensional and is spanned by the states \ipn',Q), 
{n = 1,...,A^) and IV'o;!)- Here, \tpn',0) labels a collec- 
tive atomic state in which the n^^ atom is in state all 
other atoms are in states |0), and the mode has pho- 
tons; while IV'o; 1) corresponds to all qubits in states |0) 
and one photon populating the common mode. 

To implement the required operations we exploit the 
direct coupling of the logical qubits via the common cav- 
ity photon modes, which we denote as Aj (Aj) and are 
obtained through a Fourier transform of the local modes. 
Our aim is to achieve a relatively strong coupling be- 
tween the logical qubits via the common cavity modes. 
We therefore consider a parameter regime in which the 
excited atomic states |e)fc (A; = 1, iV) evolve on a much 
faster time scale than all other states, J ^ $7, g A. 
This allows us to adiabatically eliminate these states from 
the dynamics of the system. In this case the Hamiltonian 
([1]) for qubits coupled to a single common mode A\ sim- 
plifies to 

N 

H°*' = 5]g'|0)fe(l|fc4e"''^''+H.c. + MUi|0)fc(0|fe (2) 
fc=i 

with effective coupling constants g' = —gft/2A, and de- 
tuning A' = -17V4A. Moreover, S = -g^/iA - 2J) is 
the detuning of the coupling of the qubits to the lowest 
energy Bloch mode A\. 

Performing a simple phase transformation the Hamil- 
tonian ([2]) can be rewritten in terms of the collective 



atomic states as 

N 

H<=« = ^5'|V'fc;0)(^o; i| + H.c. + (A' - J)|V'o; iXV-o; i|. 

k=l 

(3) 

Note that the above Hamiltonian describes the inter- 
action between N degenerate states, corresponding to 
\ipk',0}, coupled to one excited state \tpo'A) by coupling 
constants g' and detuned by {A' — 6). We assume that the 
qubits are addressed by a global laser pulse which auto- 
matically fulfills the requirements for equal couplings and 
detunings. For a system governed by the Hamiltonian 
from Eq. ^ Ref. [21| provides an analytical expression 
for the propagator \J{t,0) which determines the system's 
state at a moment t according to = U(t, 0)|'J'(0)). 

In the basis of the states {jV'n; 0), |?/^o; 1)} U(t,0) is given 

by 

\^{t,Q) = ^ {^,^. + (a(t)_l)^}|V;,;0)(V,;0| 

+ a*(t)|^o;i)(V'o;i| 

Here we have introduced the collective atomic Rabi fre- 
quency X — y^g'- The complex parameters a and 
b depend on the interaction between the qubits and 
the applied external fields. They can be found analyt- 
ically or numerically for any set of detunings, pulses' 
shapes and amplitudes {A' — S,g'}; and obey the rela- 
tion |ap -|- |6p = 1. Some values of interest of a and b 
are analyzed in Ref. [2l| . 

Now, we proceed with the execution of the Grover algo- 
rithm. First, we describe the initialization of the qubits 
in a \W) = 'En=i\^n)/VN = [1, . . . , l]T/y]v state, 
which corresponds to an equally weighted superposition 
of all states. It can be created in a single physical step 
just through the natural time evolution of Eq. start- 
ing from IV'o; 1) and applying a global laser pulse with a 
hyperbolic-secant time dependence, a pulse width T and 
an area of tt, while satisfying A' — 5 = 0. Hence, the 
typical width of this pulse is T = ^ = ^ , ~ 10~^s, 
for the estimation of which we have used an experi- 
mentally observed individual coupling strength g ~ 105 
MHz from Ref. j2^] and we have assumed that ft ^ g 
and A ^ lOg. For this interaction step the parameter 
a — and the state of the system thereafter is given by 

After that we continue with the description of the im- 
plementation of the Grover iteration. Mathematically, 
the Grover logical operation consists of two steps - a lo- 
cal oracle operation O which marks the searched state 
by flipping its phase and a global reflection about the 
mean. Once we have prepared the required N qubit \W) 
state the oracle operation O can be easily implemented 
by a local laser pulse, addressing only the searched qubit. 
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FIG. 2: Numerical simulation of the population of the 
searched state for the probabilistic Grover search imple- 
mented in an array consisting of A'' = 8 coupled cavities. 
The system is initially prepared in an eight-qubit W state 
through the application of a single global pulse with an area 
of TT. After the application of two Grover iterations each com- 
prising an oracle operation (implemented by a local laser pulse 
of area 2n) and a global reflection about the mean (realized 
with a global pulse of area 27r and ) the fidelity of the state 
approaches unity. After that it starts to decrease and shows 
oscillatory behaviour. The sharp disturbances of the plotted 
curve at 15T, 75T and 135r indicate the application of the 
oracle operations while the distinct increases (decreases) at 
45r, 105T and 165T show the global laser pulses which were 
chosen to have sech time dependence and a pulse area of 2n. 



with a pulse area of 27r. The resulting state of the system 
is then evolved according to the propagator given in Eq. 
(Ul as we apply a consequent global pulse, addressing all 
qubits, in order to realize the Householder reflection. It is 
again assumed to have a hyperbolic-secant shape and an 
area of 27r (hence, the individual qubits experience pulse 
areas of 2n/ ^/N) which effectively implements a 2n rota- 
tion between the states J2i=i IV'i! 0) iV'o; !)• Then, the 
parameter a = —1 and the propagator for the N degen- 
erate states (HI), comprising the database of the system, 
becomes 



g{w) ^i-2\w){w\, 



(5) 



which is the iV-dimensional Householder reflection op- 
erator about the mean [2^. Here, the normalized re- 
flection vector \W) = ^[1, 1, . . . , 1]'^ is the N qubit 

\W) state. Hence, in a single physical step we can real- 
ize in the register's subspace the basic ingredient for the 
Grover's algorithm just by tuning accordingly the inter- 
action par ameters of the system. Moreover, it has been 
shown 2JI that any desired A^-qubit unitary operation 
can be realized as a sequence of reflection operations ([5]) 
in a number of steps ~ 0{N), which is an improvement 
compared to methods based on, for example, sequences of 
two-dimensional rotations 20]. The Grover logical step 



g{W)0 is successive executed 7rviV/4 times on the reg- 
ister which drives the system into the marked state. 

The third final step in the search algorithm is the de- 
tection. The marked qubit is the only one in state |1) 
while all others are in state |0). This can be easily probed 
by employing usual atomic state measurement techniques 
|25| . In Fig. [2] a numerical calculation of the fidelity of 
the preselected state is plotted as a function of time for 
an array consisting of = 8 cavities. After the appli- 
cation of two Grover iterations, i.e. two oracle calls and 
two global reflections, almost all population is driven into 
the marked state, with fldelity 0.98%. Then, it starts to 
decrease as a part of oscillations between zero and one. 

In an experimental setup in addition to the conditions 
J ^ ri, (7 ^ A we also need to assume that the cavity 
leakage time is smaller than the time required to imple- 
ment the logical/physical steps. In our case this trans- 
lates to X < 7cav, i- e., close to strong coupling regime 
which corresponds to single-cavity cooperativity param- 

2 

eter C — — ^ 1, where T is the spontaneous emission 
rate of the atom. In addition cavities need to be effi- 
ciently coupled to each other, i.e. J > 7cav 

There are three main potential technologies for the 
experimental implementation of coupled cavity arrays: 
fiber coupled microtoroidal cavities t26|], arrays of de- 
fects in photonic crystals [131 and superconducting qubits 
coupled via microwave stripline resonators [2^. Micro- 
toroidal cavities are routinely produced in large arrays 
and can have high Q-factors [23. The cavities are cou- 
pled via tapered optical fibers which are placed close to 
the surface of the cavities and whose evanescent fields 
overlap thus allowing for photons to hop between neigh- 
bouring sites. The photon tunneling rate can be con- 
trolled by adjusting the distance between cavity and 
fiber. Moreover, these cavities can be made to inter- 
act with atoms in the strong coupling regime and single- 
cavity cooperativity parameters of C '-^ 50 have already 
been demonstrated in an experiment 30]. Some of the 
challenges for the realization of the effective Hamiltonian 
given by Eq.dH) using toroidal microcavities is that all 
the cavities of the array should be coupled and tuned 
into resonance with each other. 

Another promising candidate for implementing arrays 
of coupled cavities are atoms coupled to photonic band 
gap defect nanocavities. So far, large arrays of coupled 
cavities have been produced [3l| where photon hopping 
has been observed. Atomic impurities can be created 
inside these nanocavities and due to the cavities' small 
volume the interaction between the atom and the cavity 
is in the strong coupling regime with very large coupling 
constants [s^] and single-cavity cooperativity parameters 
up to C 10. Spontaneous emission rates can be made 
small; however, cavity decay rates of photonic band gap 
cavities remain a limiting factor for the implementation 
of the dynamics of Eq.([T|). 

Coupling of two cavity QED systems formed by a 
Cooper-pair box coupled to a superconducting stripline 
resonator has not been achieved in an experiment yet. 
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FIG. 3: Numerically calculated maximum fidelity of the 
marked state plotted against the random mean deviation of 
the different coupling constants g'i in units of 10%. 

however, these systems possess some advantages that 
make them a possible platform for the implementation 
of the Hamiltonian ((Ij. Strong coupling between the 
Cooper-pair box and the stripline resonator has been ob- 
served and very large single cavity cooperativity param- 
eters C ^ 10^ have been achieved already jSiJ]. More- 
over, two Cooper-pair boxes have been strongly coupled 
to the cavity mode [1^ HJl and theoretical investiga- 
tions suggest that their number can be increased up to 
ten. These experimental results, combined with efScient 
cavity-cavity coupling make the circuit QED system a 
promising platform for the realization of the proposed 
algorithm. 

In order to check the robustness of our protocol against 
possible experimental errors we have considered an array 
of 8 coupled microtoroidal cavities and have allowed for 
realistic inaccuracies in the resonant frequency of the cav- 
ities ujc- Another aspect which remains challenging and 
has the potential to reduce the fidelity of the protocol is 
the positioning of the atoms such that they all experi- 
ence equal coupling constants. To estimate the influence 
of these errors we have performed a numerical simulation 



of the fidelity of the marked state as a function of a static 
deviation from the desired equal value of the effective 
coupling constants experienced by the individual qubits 
g' . For the simulation shown in Fig. [3] we have assumed 
a global pulse tuned on A' — S with a hyperbolic-secant 
time dependence and a pulse area of 2tt. Our calculations 
show that the marked element is successfully recovered 
with a probability of more than 70% even when the devia- 
tions are of the order of 30%. This fidelity translates into 
success probability of about 20% for a single step, which 
is comparable with the perfect classical search algorithm 
which has a maximum average success rate of 25% for 
a single step. The error bars, however, increase rapidly 
when the mean fluctuations of the couplings are around 
20% but for smaller values they are of the order of few 
percents. This shows that we can achieve high fidelity 
with a relatively realistic system which can be crucial for 
the experimental implementation of the search algorithm. 

In conclusion we showed how symmetric cavity-qubit 
couplings can be used to implement a robust and efficient 
quantum search algorithm in a linear chain of N cou- 
pled cavities doped with atoms whose hyperfine states 
are used to encode the logical qubits. The algorithm is 
implemented by employing global laser fields in such a 
way that the number of the required physical steps be- 
come equal to the number of logical steps, i.e., ~ 0{Vn). 
This is an advantage compared to other implementations 
which in addition to these query steps require many more 
physical non-query steps. We achieve the speed up with 
regard to physical steps by using the natural evolution 
of the system, which already implements the global re- 
fiection operation required by the algorithm. We also 
allow for realistic deviations from the symmetric cavity- 
qubit couplings that can arise in an experimental setup 
and find that our proposal is quite robust against such 
imperfections. 
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